Draft version January 19, 2013 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



O 
u 

< 



X 

Oh! 

6 



> 

O 



THE COOL ACCRETION DISK IN ESO 243-49 HLX-1: FURTHER EVIDENCE OF AN INTERMEDIATE MASS 

BLACK HOLE 

Shane W. Davis^ Ramesh Narayan^, Yucong Zhu^, Didier Barret^-^, Sean A. Farrell Olivier Godet^-"', 

Mathieu Servillat^, and Natalie A. Webb'^'^ 

Draft version January 19, 2013 

ABSTRACT 

With an inferred bolometric luminosity exceeding 10^^ erg s""'^, HLX-1 in ESO 243-49 is the most 
luminous of ultraluminous X-ray sources and provides one of the strongest cases for the existence of 
intermediate mass black holes. We obtain good fits to disk-dominated observations of the source with 
BHSPEC, a fully relativistic black hole accretion disk spectral model. Due to degeneracies in the 
model arising from the lack of independent constraints on inclination and black hole spin, there is a 
factor of 100 uncertainty in the best-fit black hole mass M. Nevertheless, spectral fitting of XMM- 
Newton observations provides robust lower and upper limits with 3000 Af© ^ < 3 x lO^M©, at 90% 
confidence, placing HLX-1 firmly in the intermediate-mass regime. The lower bound on M is entirely 
determined by matching the shape and peak energy of the thermal component in the spectrum. This 
bound is consistent with (but independent of) arguments based solely on the Eddington limit. Joint 
spectral modelling of the XMM-Newton data with more luminous Swift and Chandra observations 
increases the lower bound to 6000Mq, but this tighter constraint is not independent of the Eddington 
limit. The upper bound on M is sensitive to the maximum allowed inclination i, and is reduced to 
M < 10^ Mq if we limit i < 75°. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: binaries — X-rays: 
individual (ESO 243-49 HLX-1) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are extragalac- 
tic X-ray sources whose luminosities match or exceed 
the Eddington luminosity for accretion onto a 10 Mq 
black hole (BH). It is thought that the most luminous 
objects in this class may harbor BHs with masses in the 
100 to IO^Mq range (intermediate mass black holes or 
IMBHs). BHs in this mass range are of particular in- 
terest because our current understanding of stellar evo- 
lution suggests they could not be formed by the col- 
lapse of a single ma ssive star in the current epoch of 
star formation (e.g., iBelczvnski et al.l 120101 and refer- 
ences therein). The strongest support for the IMBH 
interpretation is founded on theoretical arguments that 
BH accretion flows cannot radiate significantly above the 
Eddington luminosity. If this argument holds, the most 
luminous ULX sources, which exceed 10"*° erg s~^, must 
host IMBHs. Alternatively, if real accretion flows are 
capable of radiating significantly above the Eddington 
limit or their emission is strongly beamed ([King et al.l 
I2001f ). these sources might contain BHs of only a few 
tens of solar masses. BHs in this lower mass range have 
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been identified through dynamical studies of X-ray bina- 
ries in the Milky Way and oth er Local Group galaxies 
(|Remillard k McClintocd[2006[) . These sources can eas- 
ily be explained by the theory of stellar evolution and 
collapse. 

Cool thermal components found in the spectra of some 
ULX sources provide additional support for the IMBH 
interpretation. Due to the scaling of the characteris- 
tic gravitational radius and luminosity with mass, ac- 
cretion disks radiating at a fixed fraction of the Edding- 
ton luminosity will tend to have lower maximum effec- 
tive temperatures for larger masses if the inner radius 
of the accr etion flow scales with the gravitational radius 
of the BH (IShakura fc SunvaevllT973l : INovikov k Thornd 
[19731) . In many luminous ULXs, fits with multicolor 
disk blackbody models f avor relatively cool disks (e.g., 
IMiUer et al.|[200l [2004a llbl). However, the disk compo- 
nent in many of these spectral fits only accounts for a 
sma ll or modest fraction o f the bolometric p ower (see 
e.g. iSocrates fc Davis! I2006t iSoria et al.l 120091) . which is 
instead dominated by a power law component, generally 
thought to be inverse Compton scattering of disk pho- 
tons by hot electrons. If the majority of the emission 
originates in the hard component, it is no longer clear 
that the thermal component is associated with emission 
from near the inner edge of the disk (as opposed to larger 
radii), and the argument that a larger emitting area fol- 
lows from a larger gravitational radius (i.e. a large mass) 
is weakened. 

Making a convincing argument for an IMBH on the 
basis of spectral fits requires observations in which a 
thermal component dominates the bolometric emission. 
However, such observations appear to be rare for ULXs 
with luminosities > 10"*° erg s~^ (jSocrates fc Davi^ 
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2001: iBerehea et al l [20081: ISoria fc Kuncid [20081: 
Feng fc Kaaretll2009l: iSoria et al.ll2009D . Most luminous 
ULXs are observed in a power-law dominated state 
resembling the hard state or steep powerlaw s tate of BH 
X-ray binaries (i McChntock fc Remillarl l2006'). ahhough 
there are also suggestions that this m ay be a new mode 
of super-Eddington accreti on (e.g., iSocrates fc DavisI 
120061: [Gladstone et al.ll2009l ). E xceptions may include 
recent observations of M82 X-1 (jFeng fc Kaaretl I2010D 
and HLX-1 in ESQ 243-49, he reafter referred to simply 
as HLX-1 (jFarreh et al.ll2009[ ). in which the spectra of 
these sources appear to be dominated by a thermal 
component. 

In this work we focus on HLX-1, an off- nuclear X-ray 
source in the galaxy ESO 243-49, which reaches lumi- 
nosities in excess of 10"^^ erg s~^ jiFarrell et al.l 120091 : 
IGodet et al.[ [2001 IWebb et al.|[20Tol ). The source has 
been observed a number of times with XMM-Newton, 
Swift, and Chandra, displaying long term spectral vari- 
ability that is consis tent with state tra. nsitions in Galac- 
tic X-ray binaries (jGodet et al.l 120091 : Servillat et al., 
in preparation). This large luminosit y is con t ingent 
on its placement in ESO 243-49. So ria et al.l (|2010f ) 
have argued that the X-ray spectrum could plausibly 
be explained as a foreground neutron star in the Milky 
Way. However, based on recent spectr oscopic observa- 
tions with the Very Large Telescope, iWiersema et al.l 
([20 10;) identify Ha emission coincident with HLX-1 at 
a very similar redshift to that of the host galaxy, placing 
the source firmly in ESO 243-49. 

For fitting the spectrum of HLX-1 we use BHSPEC 
pavis et al.l 120051: [Davis fc Hubenvl I2006D . a relativis- 
tic accretion disk spectral model, w hich has been imple- 
mented as a table model in XSPEC (IArnaudlll996l ). The 
major advantage of this model is that the spectrum of the 
emission at the disk surface is not assumed to be black- 
body, but is instead c omputed directly using TLUSTY 
([Hubenv fc Lanziri995[ ). a stellar atmospheres code which 
has been modi fied to model t he vertical s tructure of ac- 
cretion disks (iHubenvl [T990t iHubenv fc Hubenv 1998). 
The BH mass (and spin) are parameters of the model 
and are directly constrained by spectral fitting. 

BH SPEC, alone or in concert with KERRBB (|Li et al.l 
|2005[) . has been used to fit numerous Galactic BH X- 
ray binaries. For observations in which the thermal 
component dominates, it provides a good fit and repro- 
duces the spectra l evolution as luminosity varies (e.g. , 
IShafee et alJ[20M [Davirerall[200l TSteiner et al.l[2010l 
and references therein). However, the mass range covered 
by the BHSPEC models used in those works was below 
the IMBH range. BHSPEC models in the IMBH regime 
were generated bylHui et all (j2005D . With these models, 
IHui fc Krolilj mm obtained good fits and found best- 
fit BH masses in the 23-73 Mq range for five of the six 
ULX sources in their sample, all of which were at least an 
order-of-magnitude lower in luminosity than HLX-1. We 
have extended our own version of the BHSPEC model 
into the IMBH regime and have used these models for 
the analysis presented he re, but the methods are identi- 
cal to those employed bv lHui fc Krol ik (2008). 

This work is organized as follows: We bricfiy describe 
the BHSPEC model and discuss our data selection in 
Section [2] We summarize the main results of our spec- 
tral analysis in Section [3l and provide a more thorough 



discussion of our constraints in Section [Jj We summarize 
and conclude in Section O 

2. METHODS 

2.1. Spectral Models 

In this work, we focus exclusively on the BHSPEC 
model. We are motivated by the observational evidence 
of a thermally dominant soft X-ray component in the 
spectra of HLX-1 considered here and the similarity of 
these spectra to those observed in Galactic BH X-ray 
binaries. Hence, we assume from the outset that the 
emission is from a radiatively efficient, thin accretion 
disk, and derive constraints on the parameters of inter- 
est, most importantly the BH mass, M. We refer readers 
interested in comparisons with other models to Godet et 
al. (in preparation), which fits a wider array of models 
and attempts to differentiate between various interpreta- 
tions. 

In addition to logM, the BHSPEC model has 4 fit 
parameters: a normalization, the BH spin a*, the cosine 
of the inclination i, and the log of the Eddington ratio 
i i/LEdd, where L^dd = 1-3 x 10^^{M/Mq) erg s'^ 
is the Eddington luminosity. The normalization can be 
fixed using the known distance to the source of D 
95 Mpc, leaving 4 parameters to fit. We have computed 
models for log^ = -1.5 to 0, log M/M© = 3.25 to 5.5, 
cosi = to 1, and a*_ = —1 to 0.99. We assume a 
iShakura fc SunyaevI ()1973n stress parameter a = 0.01, 
but discuss the implications of this choice in Section [31 

Although the spectra under consideration are ther- 
mally dominated, there is a tail of emission at high ener- 
gies which is not accounted for by BHSPEC. W e model 
this emission with SIMPL (jSteiner et al.l l2009[ ). which 
adds two free parameters: the power law index T and 
the fraction of scattered photons fsc relative to the BH- 
SPEC model. Photoelectric absorption by neutral gas 
along the line of sight is accounted for by the PHABS 
model in XSPEC. We leave the neutral hydrogen column 
Nn free, adding a single additional parameter. 

2.2. Data Selection 

We begin with a collection of prospective disk domi- 
nated observations of HLX-1. This includes two obser- 
vations with XMM-Newton (2004 November 23: XMMl; 
2008 November 28: XMM2), two observations with Swift 
(2008 November 24: SI; 2010 August 30: S2), and one 
observation with Chandra (2010 September 6). 

The XMM-Newton observations were performed with 
the three EPIC cameras in imaging mode with the thin 
filter. The two MOS cameras were in full-frame mode 
for both observations. The pn was operated in full frame 
mode for the first observation (Obs-Id: 0204540201) and 
small window mode for the second (Obs-Id: 0560180901). 
The data were reduced using the XMM-Newton Science 
Analysis System (SAS) v8.0 and event files were pro- 
cessed using the epproc and emproc tools. The event 
lists were filtered for event patterns in order to maxi- 
mize the signal-to-noise ratio against non X-ray events, 
with only calibrated patterns (i.e. single to double events 
for the pn and single to quadruple events for the MOS) 
selected. Events within a circular region of radius 60" 
around the position of HLX-1 were extracted from the 
MOS data in both observations and from the pn data in 
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TABLE 1 

BEST-FiTj^MM-Afeiuton Parameters 



cosi a, log A/'' \ogi Nh^ F fsc xVd-o.f. 

0.23t«:« 4.79t«;J0 0.001°.,, 1.51}^ 2.371^:^8 0.1081^^^ 332/328 
0.25 0.99t°.o7 5.071°;°^ -l-islo.ol 2.611:° 2.72lg:« 0.142l°:°2 339/328 
0.5 0.99l°o4 4.78lo°g -O-^Y^o'^l "^-^tl'l 2.541°;^^ 0.1081q°4j 340/328 
0.75 -0.93lj°^ 3.72lg-^® -0-16lo'°6 a.sl^j'g 2.28l°'^^ 0.080l°'°^° 338/328 

1 0.5ltlil 3.83l°-^i -0.37l°-g2 2.5li-J 2.3ll°-^9 0.085lg-°^4 335/328 
0.941°:°^ -l' 3.6ll°;i« -O.ielHi 1.9li;? 2.371°;^^ 0.0911°;°^° 336/328 
0.03l°:°o 4.70i°:?» O.OOlg.o, 2.ltli 2.38l°:« 0-099to:o22 334/328 
O.OOlo°^ 0.7 4.89lo;°7 -0-26l°:o8 ^-^-ia 2.34l°;i^| 0-099to:o26 335/328 
l.OOloig 0.95 4.07j:o;oJ, -0-43lo:i5 2.81^;^ 2.24l°;3° 0.0731q;o25 337/328 



Note. — Parameters tabulated without errors were fixed during the fit. Errors are computed 
usmg Ax^ = 2.706 (90% confidence for one parameter). Joint confidence contours are shown in 
Figures ill and |3] 

M is in units of Mq . 
^ Nil is in units of 10-^° cm~^. 




Energy (keV) 

Fig. 1. — XMM2 spectrum (top) and best-fit unabsorbed model 
(bottom) versus energy for cos j = 1 in Table [T] The top panel 
shows the observed count rate and best-fit model (with absorp- 
tion) for the EPIC pn (black), MOSl (red) and MOS2 (blue) data. 
The bottom panel shows the best-fit total model (solid, black), BH- 
SPEC component (dotted, red), and SIMPL component (dashed, 
blue) with absorption removed. 

the first observation. Background events were extracted 
for tlic same data from an annulus around tlie source po- 
sition with, an inner radius of 60" and an outer radius of 
84.85" . As the pn was in the small window mode during 
the second observation, a smaller circular extraction re- 
gion with a radius of 40" was used for extracting source 
events. Background events were in turn extracted from 
a circular region with radius 40" from a region the same 
distance from the center of the chip as HLX-1, which 
appeared to have a similar background level in the im- 
age. Source and background spectra were extracted in 
this way for each camera, with response and ancillary re- 
sponse files generated in turn using the SAS tools rmfgen 
and arfgen. 

Although the MOS and pn spectra are consistent down 
to ~ 0.5 keV, below this boundary the pn spectrum de- 
viates significantly, showing a sharp, absorption-like fea- 
ture that is absent in the MOS spectra. Since this in- 



dicates a problem with the low energy res ponse of the 
pn ca mera in this observation, we follow iFarrell et al.l 
((2009h and ignore channels with energies below 0.5 keV 
when fitting the pn spectrum. For the MOS spectra we 
fit channels with energies greater than 0.2 keV. 

The Swift XRT data were processed us ing the tool 
XRTPIPELINE vO.12.3.4, as discussed in iGodet et all 
([2009). The C/iandra ACIS-S data reduction is discussed 
in more detail in Servillat et al. (in preparation). Their 
modeling indicates that the data likely suffer fro m mild 
pile-up , so we include a pile-up model based on iDavisI 
(|2001[ ) in our spectral fits, following the guidelines in the 
Chandra ABC Guide to PileufQ- The grade migration 
parameter (denoted by a in the guide, and not to be con- 
fused with the accretion disk stress parameter) was left as 
a free parameter in our initial fits. These fits generically 
favored the maximum allowed value of 1. Our best- fit 
BHSPEC parameters were mildly sensitive to variations 
in this a, but a ^ 1 yields results which are consistent 
with the SI and S2 datasets. Hence we fixed a = 1 for 
subsequent analysis. The frame time was set to 0.8 s to 
match the observational setup, and all other parameters 
were left at their XSPEC defaults. 

All spectra were rebinned to require a minimum of 20 
counts per bin. 

Using XSPEC0 (vl2.6.0q), we fit each observation in- 
dependently to determine its suitability for modeling 
with BHSPEC. Specifically, we evaluate the level of disk 
dominance for typical best-fit parameters. For both of 
the XMM-Newton observations, we fit all EPIC datasets 
(MOSl, M0S2, and pn) simultaneously. An additional 
hard X-ray component is necessary to obtain a good fit 
for these datasets. Fits with SIMPL typically find a best- 
fit scattering fraction fsc ^ 0.5 for the XMMl data and 
fsc ^ 0. 15 for XMM2. The se results are consistent with 
those of lFarrell et al.l ()200 ty) who find an acceptable fit to 
the XMMl data with an absorbed power-law, but require 
an additional soft c omponent (modeled with DISKBB, 
iMitsuda et 311119841 ) for the XMM2 data. Our fits sug- 
gest that the XMMl spectrum is more typical of a steep 

^ http:/ /cxc. harvard.edu/ciao/download/doc/pileup_abc.pdf 
* http:/ /heasarc. nasa.gov/docs/xanadu/xspec/index. html 
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Fig. 2. — Joint confidence contours of a* and M, assuming fixed 
inclination for the XMM2 data. Each set of contours is computed 
by fitting the BHSPEC model with fixed cosi = (solid, green), 
0.25 (dotted, black), 0.5 (dashed, red), 0.75 (dot-dashed, blue), or 
1 (triple-dot-dashed, violet). The three contours in each set are 
Ax^ = 2.3, 4.61, and 9.21, corresponding to 68%, 90%, and 99% 
joint confidence on a, and logM. 

power-law state, while the XMM2 spectrum is consis- 
tent with a thermally dominated state. A large fsc in 
the fit with SIMPL is problematic for our accretion disk 
model, which does not account for the effects of irradi- 
ation by neighboring corona. Hence, we only report our 
results from the XMM2 analysis below. We note that the 
XMMl data generally favor a best-fit M that is higher 
than XMM2 for the same i. This yields a cooler disk 
model, consistent with the fact that more of the high 
energy flux is accounted for with the SIMPL component 
in these data. However, the joint confidence contours on 
a* and M are much larger for the XMMl data, so the 
best-fit values are still consistent with the XMM2 values 
at 90% confidence. 

The Swift and Chandra observations caught the source 
at higher luminosities than the XMM2 observation, but 
the overall signal-to-noise is still lower because of the 
lower effective area of the Swift XRT and the short dura- 
tion (10 ks) of the Chandra ACIS-S observation. Due to 
the lower signal-to-noise and stronger disk dominance, a 
suitable fit is provided by the (absorbed) BHSPEC model 
alone. The addition of the SIMPL model only provides a 
marginal improvement to the fit and the best-fit SIMPL 
parameters are poorly constrained. We also fit the com- 
bined SI and S2 observations. In this case we tie all 
BHSPEC parameters together, except for £, which is al- 
lowed to vary independently for SI and S2. Even for 
the combined dataset, we find a good fit with BHSPEC 
alone, and poor constraints on the SIMPL model param- 
eters. Therefore, we do not include the SIMPL model in 
subsequent analysis of the combined SI and S2 datasets 
or in our analysis of the Chandra data. 

3. RESULTS 

3.1. XMM-Newton Results 

We first consider the XMM2 observation. The soft 
thermal component is largely characterized by only two 
parameters: the energy at which the spectrum peaks and 
the overall normalization. In practice, this leads to de- 
generacies in the best-fit parameters of BHSPEC (see 



Fig. 3. — Joint confidence contours for BH mass and inclination, 
assuming fixed BH spin, for the XMM2 observation. Each set of 
contours is computed by fitting the BHSPEC model with fixed 
a» = —1 (solid, green), (dotted, black), 0.7 (dashed, red), or 0.95 
(dot-dashed, blue). The three contours in each set correspond to 
68%, 90%, and 99% joint confidence on cosi and log A/. 

e.g. lDavis et al.ir2006f) . unless all but two of the parame- 
ters can be independently constrained. Since we can only 
constrain D independently, we will need to consider joint 
variation of the remaining parameters. 

For completeness, our best-fit parameters for various 
choices of i or a* are summarized in Table [T] We find 
acceptable fits for all i. The best- fit model and data 
for i = are plotted in Figure [TJ We report the 90% 
confidence interval for a single parameter, but due to 
the significant degeneracies in the model, the joint con- 
fidence contours better illustrate the actual parameter 
uncertainties. Hence they are the focus of this work. 

Since M and a* are of primary interest to us, we first 
examine their joint confidence contours, and consider the 
joint confidence of M and i in SectionlH For illustration, 
it is useful to consider several sets of contours for different 
choices of (fixed) i, leaving £ as a free parameter. These 
best-fit joint confidence contours are shown in Figure [51 
We consider five choices for cosi, evenly spaced from to 
1. Each set of contours has three curves corresponding 
to 68%, 90%, and 99% confidence. These are determined 
by the change in relative to the best-fit values listed 
in Table □ 

Even for fixed i, there is a clear correlation of a* 
and M. At 99% confidence, the entire range of a* 
(— 1 < a, < 0.99) is allowed, and the corresponding 
M varies by a factor of 4-5 over this range. The only 
exception is for nearly edge on systems (cosi 0), for 
which low spins are disallowed. The correlation of best- 
fit M with i is very strong as well, c onsistent with previ- 
ous modeling of other ULX sources (IHui fc Krolikl[20081 
This is illustrated more clearly in Figure |31 which shows 
the joint confidence contours of M and cosi for fixed 
a* = —1, 0, 0.7, and 0.95. The corresponding best-fit 
values are listed in Table [TJ At fixed a* , there is up to 
a factor of 10 change in AI as inclination varies from 0° 
to 90°. The combined overall uncertainty in the mass is 
nearly a factor of 100. 

Despite these degeneracies, we can still place a firm 
lower limit of M > SOOOM© on the mass of the BH in 
HLX-1. The limiting mass is obtained for the case of 
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Fig. 4. — Prom top to bottom, the best fit Eddington ratio, 
power law index, and scattering fraction in the a, — log M plane, 
assuming cosi = 0.5. For comparison, we overplot three contours 
corresponding to 68%, 90%, and 99% joint confidence on a* and 
logM. 

a maximally spinning BH with a counter-rotating disk 
viewed nearly face-on. An independent lower limit on 
the mass may be obtained under the assumption that the 
source luminosity does not exceed the Eddington limit. 
This gives the same answer, M > SOOOM©, which is 
purely coincidental. At the high mass end, the limit we 
obtain from our fits is M < 3 x 10^ Mq, where the lim- 
iting value is obtained for nearly maximal spins in the 
case where disk rotation is spin-aligned and the system is 
viewed nearly edge on. However, edge on systems may be 
ruled out by the lack of X-ray eclipses if one assumes the 
accreting matter is provided by a binary companion. If 
one limits i < 75°, one obtains M < 10^ M©. The above 
limits on the mass place HLX-1 in the IMBH regime, 
though the very highest masses in our allowed range ap- 
proach the lower end of the mass distribution inferred in 
active galactic nuclei (e.g., Greene & Ho 2007). 

We plot the variation of the other free parameters in 
Figure m for the specific example of cosi = 0.5. The 
variation of £ is the most interesting, as there are fairly 
clear correlations with M and a*. Following the best-fit 
contour, £ decreases as M and a* increase. This cor- 
relation arises because £ has a role both in determining 
the location of the spectral peak as well as its overall 
normalization, a point we elaborate on in Sectional 

The bottom two panels show the variation of the 
SIMPL parameters. In contrast to £, these parame- 
ters tend to vary across (rather than along) the confi- 
dence contours, indicating that they are highly corre- 
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Fig. 5. — Joint confidence contours of a* and At, assuming fixed 
incfination for the Swift data. Each set of contours is computed 
by fitting the BHSPEC model with fixed cosi = (solid, green), 
0.25 (dotted, black), 0.5 (dashed, red), 0.75 (dot-dashed, blue), or 
1 (triple-dot-dashed, violet). The three contours in each set are 
Ax^ = 2.3, 4.61, and 9.21, corresponding to 68%, 90%, and 99% 
joint confidence on a* and logAf. 

lated with x^. A similar behavior is present in the ab- 
sorption column. For the best-fit values we find Nn ~ 
1 - 4 X 10^° cm^2, r ~ 2.2 - 2.7, and a scattered fraction 
< 12%. The allowed iVn range is roughly equivalent to or 
slightly gre ater than the Galactic value (1.7 x 10^° cm~^, 
IKalberla et al. 2005). The ranges for the two SIMPL pa- 
rameters are generally consistent with fits to the thermal 
state of Galactic BH binaries and confirm that the bolo- 
metric luminosity of the best-fit models is dominated by 
the BHSPEC component. For typical parameters within 
the confidence contours, the accretion disk accounts for 
80-95% of the unabsorbed model luminosity from 0.3 to 
10 keV . Comparison with Table 1 of Socrates & DayJi 
(|2006f l shows that this is a much higher disk fraction 
than typically inferred in other bright ULXs, in which 
values 10-50% are more typical. 

3.2. Swift and Chandra Results 

In Figure [SJ we plot the best-fit joint confidence con- 
tours of a* and M for the combined Swift (SI and S2) 
data sets. Figure [6] is the corresponding plot for the 
Chandra dataset. As discussed in Section 12.21 we do 
not include a SIMPL component in fitting these datasets 
and we tie all the BHSPEC parameters together, except 
£, which is allowed to vary independently for each ob- 
servation. The corresponding best-fit parameters are de- 
picted with crosses and summarized in Tables [2] and [3l 
where £i and £2 correspond to the SI and S2 observa- 
tions, respectively. The reduced for the Chandra data 
indicate a slightly poorer fit than with the XMM-Newton 
or Swift data. This is primarily due to narrow residuals 
near 0.6, 1.1 and 1.3 keV in fits to the ACIS-S data. Since 
these residuals are not highly significant and their pres- 
ence does not significantly impact the best-fit BHSPEC 
parameters, we do not attempt to model them with ad- 
ditional components. Further discussion can be found in 
Servillat et al. (in preparation). 

The Swift and Chandra data seem to place much 
tighter overall constraints on a* and M, with a mini- 
mum a* at each i. The overall minima of M ~ GOOOM© 
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TABLE 2 
Best-Fit Swift Parameters 



cos i 


a* 


log A/'' 


log^i 


log ^2 




xVd.o.f. 




0.25 
0.5 
0.75 
1 


r, OI +0.18 

'^■'^^-0.02 

46+° *^ 
^■^°-0.07 

u. J»_Q 25 

0-99^0.05 


4 96+0.30 
^•^"-0.08 
4 4Q+0.38 
^•^^-0.08 
4 21+0.63 

4 54+0.09 
-0.58 


-0 01+" °^ 

-0 02+°-°^ 
'^•'-'^-0.54 

-0 29+"-^^ 
^•^^-0.04 


-0 05+0.03 

-0 06+° °=' 
'-'■""-0.56 

-0 06+° °"' 
-0 33+0.31 
p, ns+o.os 

-O.US_o 06 


n 9+3.1 
'-'•^-2.6 

5 3+3.5 

c: 0+3.5 
o.o_2 y 


58.3/75 
59.6/75 
59.8/75 
58.9/75 
58.5/75 



Note. — Parameters tabulated without errors were fixed during the fit. 
Errors are computed using Ax^ = 2.706 (90% confidence for one parameter). 
Joint confidence contours are shown in Figure [5] In these joint fits to the SI 
and S2 data sets all BHSPEC parameters are fit simultaneously to both spectra, 
except £ which is fit independently. Hence, £i and £2 correspond to the SI and 
S2 observations, respectively. 



^ M is in units of Mq . 

^ Nil is in units of 10^° cm.-^. 

for Swift) and and M ^ 4000Mo for Chandra both cor- 
respond to i = 0°. Fits with retrograde accretion flows 
(a* < 0) are inconsistent with the Swift data and models 
with a* < —0.5 are ruled out by the Chandra spectrum. 
For values of a* allowed by both data sets, the position 
of the confidence contours for a given i are consistent 
with each other. They are also broadly consistent with 
the XMM2 results, ahhough offset to slightly higher M. 
However, due to the lower signal-to-noise of the Swift and 
Chandra datasets, the range of allowed M for a given i 
is larger than in Figure [21 

The differences between the XMM-Newton fits and the 
Swift or Chandra results at low (negative) a, are driven 
primarily by the larger luminosities of the disk compo- 
nent in the Swift and Chandra observations. With the 
exception of nearly edge on systems (cosi = 0) we find 
£ < 1 within the 99% confidence contours in Figure [2l 
In the Swift and Chandra data, £ > 1 is required for 
all i, at sufficiently low a». As we discuss further in 
Section m reductions in a, and M must be offset by in- 
creases in £, but our models are capped at £ — 1, which 
sets a lower limit on the allowed a* and M for a given 
choice of i. Hence, in contrast to the XMM2 observation, 
these tighter constraints depend strongly on the assump- 
tion that luminosity does not exceed the Eddington limit: 
£<1. 

We also fit the combined Swift, Chandra and XMM- 
Newton datasets assuming no correction was needed to 
the relative effective area of the various instruments. In 
principle, a constant offset coul d be included based o n 
cross-calibration analysis (e.g., iTsuiimoto et al.l I2011D . 
but even without such a correction, we obtain good fits to 
the combined dataset. The range of allowed a* is similar 
to that found with the Swift data alone, but the width 
of the contours are comparable to those from the XMM2 
fits plotted in Figure [21 At 99% confidence, the mini- 
mum allowed M was about 6OOOM0, again for i = 0° 
and £ = 1. 

4. DISCUSSION 

Due to the complex interplay between the frequency 
and angular distribution of emitted photons in the rest 
frame of the flow, and the relativistic effects on the disk 
structure and radiation, it is preferable to use a self- 
consistent relativistic model to draw quantitative conclu- 
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= 0.75 
= 1.00 
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4.5 

log (M/M„ 



5.0 



5.5 



Fig. 6. — Joint confidence contours of a* and M, assuming fixed 
inclination for the Chandra data. Each set of contours is computed 
by fitting the BHSPEC model with fixed cosi = (solid, green), 
0.25 (dotted, black), 0.5 (dashed, red), 0.75 (dot-dashed, blue), or 
1 (triple-dot-dashed, violet). The three contours in each set are 
Ax^ = 2.3, 4.61, and 9.21, corresponding to 68%, 90%, and 99% 
joint confidence on a* and logM. 

sions. Nevertheless, one can understand the basic trends 
and correlations with a simpler model, in which the rel- 
ativistic effects are encapsulated into two parameters, 
similar in spirit to the pioneering work of iZhang et alj 
(fT997l) . 

4.1. A Minimal Relativistic Disk Model 

To first approximation, thermal state disk spectra are 
characterized by only two parameters: a peak energy and 
an ove rall normalization. I ndeed the popular DISKBB 
model (jMitsuda et al.l[T98^ . with only two parameters, 
provides an acceptab le fit to the soft the rmal component 
of the XMM2 data ([Farrell et al.ll2009[ ). Although we 
obtain all our constraints by fitting a relativistic accre- 
tion disk model directly to the spectrum, our fit results 
are easily understood in terms of these two constraints. 
These observables can be mapped onto a characteristic 
color temperature Tobs (~ from DISKBB), and since 
the distance D is known, a characteristic luminosity Lobs- 

The characteristic temperature is related to the peak 
effective temperature T^s of the accretion disk. For il- 
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TABLE 3 
Best-Fit Chandra Parameters 



cos i 




log M'^ 


log£ 








0.25 
0.5 
0.75 
1 


7n+0-2T 

"^-0.05 
11 +0-85 

0.701?:^^ 

0.831°;}^ 


A qn+0-28 

*-yu_o.o4 

4 07+0.37 
4 0fi+°-" 
^•^^"'-0.39 


-0 m+o-"! 
-0 ni+o oi 
00+"°° 

-CI 25+°-25 

-0 ifi+o if* 


n o+^-t^ 
1 q+3.5 

^^■•^-1.3 
1 S+4-0 
^■«-1.8 
9 0+4.8 
^■•^-2.3 
9 (1+5.8 
^■'^-2.0 


48.4/38 
50.4/38 
50.1/38 
49.4/38 
49.9/38 



Note. — Parameters tabulated without errors were fixed during 
the fit. Errors are computed using Ax^ = 2.706 (90% confidence for 
one parameter). Joint confidence contours arc shown in Figure |6] 

M is in units of Mq. 
^ TVh is in units of lO'^'' cm~'^ . 

lustration, we assume this is equivalent to the effective 
temperature at the inner edge of the disk. Ignoring rela- 
tivistic terms and constants of order unity, the flux near 
the inner edge of an accretion disk can be approximated 

by 

GMM 



(1) 



where G is Newton's constant, and M is the accretion 
rate. The inner radius of the disk corresponds to the 
innermost stable circular orbit (ISCO) in our model. 

Due to deviations from blackbody emission and rela- 
tivistic effects on photon propagation, Tobs 7^ ?eff- We 
parametrize the deviations from blackbody via a spectral 
hardening factor (or color correction) /coi and the rela- 
tivistic energy shifts with S so that Tobs = fcoiSTgg. In 
practice, 6 is primarily a function of i and a, with 5 > 1 
typical for most a* and i in our models (i.e. Doppler 
blueshifts are generally more important than Doppler 
and general relativistic redshifts). In contrast, /coi is 
dependent on M, a*, and £ and largely independent of i. 
Using this parametrization and equation ([T]) we obtain 



Tobs — To /col (5 



1/4 



— Rin/R: 



9' 



(2) 



where, m 

(c^/GA/qKosO")^^^, Rg — GM/c^ is the gravitational ra- 
dius, and Kos is the electron scattering opacity. Here, 
(. = TjMc^ / L^dd, and for simplicity we have approxi- 
mated 7] ~ 1 /"Tin • 

Using the above definitions the observed luminosity is 



Lobs — Lolmfi, 



(3) 



where Lq = AttGMqc/ k^s is the Eddington luminosity 
for a one solar mass BH and /i is a variable that en- 
capsulates all of the angular dependence of the radia- 
tion field. In an isotropically emitting, Newtonian disk 
fi = cosz, accounting for the inclination dependence of 
the disk projected area. In our models, two other ef- 
fects are important as well: electron scattering and rela- 
tivistic beaming which tend to make the disk emission 
more limb-darkened and limb-brightened, respectively. 
The latter depends significantly on the spin so fi can 
be a strong function of well as i. Lensing by the 

BH can also be important, particularly for nearly edge 



cs 0.0 



2.2 



1.8 



1.4 



1.0 

3.6 3.8 4.0 4.2 4.4 4.6 4.8 5.0 

log (M/M,„J 

Fig. 7. — Color correction corresponding to the best-fit model 
in the a, — log A/ plane, assuming cosi = 0.5. For comparison, we 
overplot three contours corresponding to 68%, 90%, and 99% joint 
confidence on a* and logM. 



on systems where it places a lower limit on the effective 
projected area of the disk. 

Equations 1^ and ^ can be solved for m and £, yield- 
ing 

/ rpi T \ 1/2 1-2 r2 

M^Mo(i^) (4) 




The quantities in parentheses on the right hand side of 
equation (j4]) are observational constraints or constants 
of nature. Uncertainties in D enter through Lobs oc D^, 
so we approximately have M cx D~^. Therefore, the rel- 
atively small uncertainty in D contributes only a very 
modest additional uncertainty in M. The remaining 
quantities arc functions of model parameters: rin(a*), 
^(a*,i), S{a^,i), and fco\{a.^,£,M). 

4.2. Correlations among the Best-fit Parameters 

The dependence of rin on a* is the primary driver of 
the correlation of M with a* in Figure [S] As a* increases 
from -1 to 0.99, nn decreases from 9 to 1.45, driving cor- 
responding increases in M . This variation is strongest 
as a* — )■ 1, leading to the "bend" in the confidence con- 
tours at high a*. Both 5 and /coi increase with a^,, and 
contribute to the correlation as well. 

The spectrum in BHSPEC is calculated directly, so 
there is no explicit /coi, but we can estimate /coi by tak- 
ing the best fit BHSPEC models and fitting them with 
the KERRBB model, for which /coi is a parameter0 The 
variation of /coi is shown for cosz = 0.5 in Figure [T] Al- 
though /col can vary by more than a factor 2 over the 
parameter range of interest, its change within the best 
fit contours is significantly less. For cosi ~ 0.5, 6 in- 
creases modestly with a* due to Doppler blueshifting, so 
the product of /coi and Toff must decrease to maintain 
agreement with Tobs- Since Toff is the primary driver of 
variations in /coi, and fcoi and T^s are positively corre- 
lated, both Toff and /coi must decrease as a^, increases. 

^ Specifically, we fix a,, i, M to correspond to the BHSPEC 
parameters, leaving /coi and M as free parameters in the KERRBB 
fit. 
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The same argument applies outside the best-fit contours, 
but leads to a much larger variation in /coi- 

The confidence contours in Figure [3] confirm a strong 
anti-correlation of M with cosi, or equivalently, a corre- 
lation of M with i. This correlation is driven primarily 
through the i dependence of 6 and fi. As i increases, the 
projected area of the disk decreases while limb darken- 
ing shifts intensity to lower i. The combined effects lead 
to a significant decrease in fi as i increases. For higher 
a*, these effects are somewhat mitigated by relativistic 
beaming, which tends to shift intensity to larger i. In ad- 
dition, Doppler blue-shifts due to the Keplerian motion 
cause S to increase with i. Both effects are present and 
drive a positive correlation of M and i, consistent with 
Figure |3l The anti-correlation of i and fi dominates at 
low spin, while the correlation of 6 with i dominates at 
higher spin because rotational velocities are a large frac- 
tion of the speed of light. The effects contribute compa- 
rably for a* ^ 0.7. 

We note that iHui fc KrolikI ()2008[ ) found a similar cor- 
relation in their analysis of several ULX sources. To 
explain their results, they attribute the correlation to 
Doppler shifts (i.e. the dependence of 6 on i), which is 
consistent with the fact that their best-fit models favored 
high a*. 

The anti-correlation of best-fit £ with M and a* follows 
directly from equations ([2]) and ([3]). As M increases, £ 
must decrease to keep Lobs approximately constant, but 
keeping Tobs constant then requires an increase in a* (i.e. 
a reduction in rin). This is why high M and low a* (and 
vice- versa) yield poor fits for cosi — 0.5 in Figure SI 

With this understanding of how M, i, £, and a* cor- 
relate, it is useful to consider what ultimately sets the 
minimum allowed M in the two data sets. The Lobs con- 
straint (eq. [3]) allows M to decrease as long as there is 
a corresponding increase in /i or an increase in i. For fits 
to the XMM2 data, i is set by the Tobs constraint (eq. 
[5]). Since £ increases as M decreases, rin must increase. 
Hence, the maximum £ and minimum M are obtained 
for a* = —1. This argument leads to a minimum M for 
any fixed i and since fi is maximum for face on disks, the 
'global' minimum of M occurs at i = 0° and corresponds 
to £ slightly less then unity. Hence, the constraints that 
M > 3OOOM0 is independent of the argument that £ < 1. 

For the Swift and Chandra datasets, which consist of 
observations with higher iobs, this is not the case. Our 
models are capped at £ — 1 due to inconsistencies in the 
underlying model assumptions for £ ^ I. For some min- 
imum M {M - 6000Mq: Swift; M - 4OOOM0: Chan- 
dra), equation ^ cannot be satisfied for £ < 1, even with 
1 = 0°. This corresponds also to a lower limit on a* since 
further increases in rjn cannot be offset by decreases (in- 
creases) in M {£) to keep Tobs constant. If we ignored 
the internal inconsistencies in thin disk assumptions and 
extended our models to £ > 1, we expect we would ob- 
tain reasonable fits for lower M and a* as equations ^ 
and ([3]) could still be satisfied. Hence, the more stringent 
lower limit on M and a* for the Swift and Chandra fits 
is not independent of the Eddington limit. 

For the XMM2 data, the minimum allowed M is 
reached for models with £ ~ 0.7, just below the Ed- 
dington limitH The assumptions underlying the thin 

^'^ Our £ only accounts for the luminosity of the BHSPEC com- 



accretion disk model are likely to break down for £ ~ 1 
although it is difficult to precisely estimate the £ value 
where our spectral constraints are no longer reason- 
able. Gene ral relativistic magnetohydrodynamic sim ula- 
tions (e.g., iPenna et al.' 2010; Kulka rni et al.l 120111 ) and 
slim disk models (ISad owski et al. 2010() suggest that the 
thin disk spectral models remain reliable to (at least) 
£ < 0.3. Observations of Galactic X-ray bi naries suggest 
that there are no abrupt changes at this £ (jSteiner et al.l 
IMS), so it is plausible that models remain basically 
sound even for somewhat higher £. 

The strong correlation of with F and fsc in Figure |4] 
(as well as a similar one for iVu) results from the failure 
of the BHSPEC model to adequately approximate the 
thermal emission for M and a* outside the confidence 
contours. If BHSPEC is a poor match to the thermal 
spectrum, the fit compensates by adjusting A^h or the 
SIMPL parameters. For example, at low M and a,, the 
disk is too cold to match the XMM2 data, so SIMPL 
adjusts by increasing the scattering fraction and making 
the power law steeper, to better fit the high energy tail 
of the thermal emission. To prevent an excess at lower 
energies, the A^h increases simultaneously. For the lower 
signal-to-noise Swift and Chandra data, A'h adjusts in a 
similar manner, even though the SIMPL component is 
absent. The width of the confidence contours appears 
to be set largely by the effectiveness of these compen- 
sation mechanisms. Hence, the extent of the confidence 
contours for a given i could presumably be reduced with 
better statistics at high energies to constrain the SIMPL 
parameters and independent constraints on A^h- 

Finally, it is worth considering the degree to which 
our correlations could be reproduced by a non-relativistic 
analysis, but assuming rjn is equal to the IS CO radius. 
For example, assuming constant /coi ~ 1.7, 6 ^ 1, and 
/i cos i in equation |4j one can recover some aspects 
of the inferred correlations of M with a* and i. The 
sensitivity of M to a, would be well approximated for 
low-to- moderate a», but the dependence of S, /coi, and 
on a, can lead to modest discrepancies as a» — 1. In con- 
trast, the correlation of M with i would only be crudely 
reproduced since the above prescription would suggest 
M oc (cosi)~^/^. This underestimates the sensitivity of 
M to i for low to moderate i, but overestimates it as 
i — >■ 90°, where the projected area goes to zero in a non- 
relativistic model. Since the low to moderate range is 
probably more relevant for observed systems, the overall 
uncertainty of M would be underestimated. We also note 
that although the assumption of a constant /coi ~ 1.7 
turns out to be approximately correct (/coi ~ 1.8 ± 0.1 
in Figure [7]), it was not clearly justified a priori and may 
not be as good of an assumption for other ULX sources. 
In principle, one could improve this analysis by estimat- 
ing /col, and 6 from BHSPEC (or some similar model), 
but at that level of sophistication, it seems more sensible 
to fit the relativistic model directly. 

4.3. Sensitivity to BHSPEC Model Assumptions 

ponent, the total bolometric luminosity Lbol of the systems is 
larger once the coronal component is added. For typical SIMPL 
parameters, the coronal component accounts for < 25% of the to- 
tal model flux, even after extrapolation to larger energies. Hence, 
ibol/i'Edd ^ 1 still holds. 
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There are a number of assumptions present in the BH- 
SPEC model that could have some impact on /coi- In 
particular, magnetic fields (and associated turbulence) 
may play a role in modifying the d isk vertical struc- 
ture an d radiative trans fer (e.g. se e Davis ct al. 2005'; 
Davis fc Hubenv ,200a IBlaes et all ^200.6.; Davis c t al.| 
2009() . Another assumption of interest is our choice 
of a = 0.01. For a < 0.01, the models depend very 
weakly on a for the parameter range relevant to our fit 
results. For higher values of a, the typical color correc- 
tion is larger. For low to moderate £ and a,, /coi in- 
creases by less than 25 % as a increases from 0.01 to 0.1. 
Much larger shifts can occur if both £ and a* are larger 
(a, > 0.8 and i > 0.3; iDone fc Dav"ii [20081 ) . but models 
in this range overpredict Tobs and are irrelevant to our re- 
sults. Hence, if the characteristic a associated with real 
accretion fiows is larger (as some models of dw arf no- 
vae a nd some numerical simulations suggest, iKing et al.l 
|2007[) . the effect would be to shift our best-fit contours 
to higher M, but only by a modest amount. For the 
parameters corresponding to the lower M limit {i = 0°, 
€ = 0.7, and a* = -1), BHSPEC yields /coi - 2. From 
equation (jU we see that reducing /coi — 1 (the absolute 
minimum) only reduces M by a factor of 4, still placing 
HLX-1 in the IMBH regime. 

Alternatively, one could make the disk around a low M 
BH look cooler by truncating it at larger radius. Equa- 
tion (U) suggests that decreasing M to a value near 30Mq 
would require a factor of 100 increase to nn ^ 900. Such 
an interpretation would need to explain why the flow 
does not radiate inside this radius. Since the energy does 
not come out in the hard X-rays, it cannot be a transition 
to an advection dominated accretion flow, which is often 
invoked to explain the low state of Galactic X-ray bina- 
ries (e.g.,.Esin et al.!!2001). Furthermore, since 77 ^ l/^jn, 

the required M would increase by a factor of 100 to 
M - 10-2 Mq yr-i. For M ~ lO"'' Mq yr-\ the ac- 
cretion rate and time variability of HLX-1 p resent a chal- 
lenge to standard models of mass transfer (jLasota et al.l 
|2011|) . Hence, it is unlikely that such a high rate is even 
feasible in a binary mass transfer scenario. 

Finally, one could plausibly obey the Eddingtoii limit 
by assuming a large b eamin g facto r Ce.g. IKing et al.|[200ll : 
IFreeland et ani2006t IKind l2008D so that iobs > iiso, 
which (in this formalism) is equivalent to increasing /i 
for some narrow range of i. Obeying the Eddington limit 
with M ~ 3OM0 would require /i ^ 100, but note that 
this is insufficient for explaining Tobs due to the /i^^^ 
dependence in equation Fixing £ = I and decreas- 
ing M by a factor of 100 yields a factor of 3 increase 
in Tobs- Maintaining agreement with Tobs requires M 
and £ to decrease proportionately, which requires a fac- 
tor of 1000 increase in fi. Any scenario with such large 
beaming factors probably requires a relativistic outflow 
or very different accretion flow geometry so these simple 
scalings may not strictly apply. Nevertheless, we empha- 
size that explaining the soft emission in HLX-1 presents a 
serious challenge to any beaming model, but is naturally 
explained by the IMBH interpretation. 

5. SUMMARY AND CONCLUSIONS 

Using BHSPEC, a fuUy-relativistic accretion disk 
model, we fit several disk dominated observations of 



HLX-1 for which the luminosity exceeds 10 erg s~^. 
Due to degeneracies in the best-fit model parameters, 
90% joint confidence uncertainties are rather large, yield- 
ing a factor of 100 uncertainty in the best-fit BH mass. 
For fits to the XMM-Newton data, we obtain a lower 
limit of M > 3OOOM0, where the limit corresponds to 
i — 0°, a, = —1, and £ — 0.7. We emphasize that 
this limit is driven by the need to reproduce the shape 
and peak energy of the thermal component in the spec- 
trum. Hence, the Eddington limit plays no role in this 
constraint. Constraints from fits to Swift and Chan- 
dra observations, which correspond to higher luminosi- 
ties, nominally offer a more restrictive lower bound of 
M > 6000 Af©, but this bound is subject to the Edding- 
ton limit because our model grid is limited to a maximum 
luminosity £ = 1. 

We also find an absolute upper bound of A/ < 3 x 
IO^Mq with both datasets, this limit corresponding to 
nearly edge-on {i — 90°) disks with near maximal spins 
(a* ~ 0.99). This upper limit is subject to the uncer- 
tainties in the models at very high spin and high in- 
clination (most notably our neglect of returning radi- 
ation and assumption of a razor thin geometry). The 
lack of X-ray eclipses and the absence of evidence for 
nearly edge-on X-ray binary s ystems in the Milky Way 
(|Naravan fc McClintocU [20051) motivate a limit on i < 
75° and, therefore, M < 10^. An argument against 
i ~ 90° based on absence of eclipses assumes that the 
accreting matter is being provided by a binary compan- 
ion, but obscuration by a flared outer disk may generi- 
cally limit the range of observable i. For M > 10^ Mq, 
HLX-1 would be consistent with the lower end of the 
mass distribution in ferred in active galactic nuclei (e.g., 
iGreene fc Hoi [200 7!), but would still be distinctive be- 
cause of its off-nuclear location in ESO 243-49. 

Other parameters of interest, such as i and a*, are es- 
sentially unconstrained by the data, unless we require 
that the disk must radiate below the Eddington lumi- 
nosity, in which case a* > or a* > —0.5 are required 
by the Swift and Chandra data, respectively. Observa- 
tions with improved signal-to-noise are unlikely to sig- 
niflcantly tighten these M constraints, as the allowed M 
range is set primarily by uncertainties in i and a* . Inde- 
pendent estimates for a* and i are ultimately needed to 
improve our M constraints, and could plausibly be pro- 
vided b y modeling o f broad Fe Ka lines or X-ray polar- 
ization ()Li et al.l 12009) if such data became available. If 
a broad Fe line is present in HLX-1, obtaining the signal- 
to-noise necessary to resolve it would require unfeasibly 
long exposure times with XMM-Newton and other exist- 
ing X-ray missions. However, such constraints may be 
possible for future missions with larger collecting areas. 

In summary, despite the rather large range of M al- 
lowed by the spectral fits of HLX-1 presented here, our 
study strongly suggests that the BH in HLX-1 is a gen- 
uine IMBH. 



We thank the anonymous referee for useful comments 
that improved the paper. SWD is supported in part 
through NSERC of Canada. SAF gratefully acknowl- 
edges funding from the UK Science and Technology Re- 
search Council and the Australian Research Council. RN 
acknowledges support from NASA grant NNXIIAEI6G. 



10 



MS acknowledges supports from NASA/Chandra grant 
DD0-11050X and NSF grant AST-0909073. This re- 
search has made use of data obtained from the Chandra 



Data Archive and software provided by the Chandra X- 
ray Center. Based on observations from XMM-Newton, 
an ESA science mission with instruments and contribu- 
tions directly funded by ESA Member States and NASA. 



REFERENCES 



Arnaud, K. A. 1996, in Astronomical Society of the Pacific 
Conference Series, Vol. 101, Astronomical Data Analysis 
Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17 — h 

Belczynski, K., Bulik, T., Fryer, C. L., Ruiter, A., Valsecchi, F., 
Vink, J. S., & Hurley, J. R. 2010, ApJ, 714, 1217 

Berghea, C. T., Weaver, K. A., Colbert, E. J. M., & Roberts, 
T. P. 2008, ApJ, 687, 471 

Blaes, O. M., Davis, S. W., Hirose, S., Krolik, J. H., & Stone, 
J. M. 2006, ApJ, 645, 1402 

Davis, J. E. 2001, ApJ, 562, 575 

Davis, S. W., Blaes, O. M., Hirose, S., & Krolik, J. H. 2009, ApJ, 
703, 569 

Davis, S. W., Blaes, O. M., Hubeny, I., & Turner, N. J. 2005, 
ApJ, 621, 372 

Davis, S. W., Done, C, & Blaes, O. M. 2006, ApJ, 647, 525 

Davis, S. W., & Hubeny, I. 2006, ApJS, 164, 530 

Done, C, & Davis, S. W. 2008, ApJ, 683, 389 

Esin, A. A., McClintock, J. E., Drake, J. J., Garcia, M. R., 

Haswell, C. A., Hynes, R. I., & Muno, M. P. 2001, ApJ, 555, 

483 

Farrell, S. A., Webb, N. A., Barret, D., Godet, O., & Rodrigues, 

J. M. 2009, Nature, 460, 73 
Feng, H., & Kaaret, P. 2009, ApJ, 696, 1712 
— . 2010, ApJ, 712, L169 

Freeland, M., Kuncic, Z., Soria, R., & Bicknell, G. V. 2006, 

MNRAS, 372, 630 
Gladstone, J. C, Roberts, T. P., & Done, C. 2009, MNRAS, 397, 

1836 

Godet, O., Barret, D., Webb, N. A., Farrell, S. A., & Gehrels, N. 

2009, ApJ, 705, L109 
Greene, J. E., & Ho, L. C. 2007, ApJ, 670, 92 
Hubeny, I. 1990, ApJ, 351, 632 
Hubeny, I., &c Hubeny, V. 1998, ApJ, 505, 558 
Hubeny, I., k Lanz, T. 1995, ApJ, 439, 875 
Hui, Y., k Krolik, J. H. 2008, ApJ, 679, 1405 
Hui, Y., Krolik, J. H., & Hubeny, I. 2005, ApJ, 625, 913 
Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., 

Bajaja, E., Morras, R., & Poeppel, W. G. L. 2005, VizieR 

Online Data Catalog, 8076, 
King, A. R. 2008, MNRAS, 385, L113 

King, A. R., Davies, M. B., Ward, M. J., Fabbiano, G., & Elvis, 

M. 2001, ApJ, 552, L109 
King, A. R., Pringle, J. E., & Livio, M. 2007, MNRAS, 376, 1740 



Kulkarni, A. K., et al. 2011, ArXiv e-prints 

Lasota, J. -., Alexander, T., Dubus, G., Barret, D., Farrell, S. A., 

Gehrels, N., Godet, O., & Webb, N. A. 2011, arXiv:1102.4336 
Li, L.-X., Narayan, R., & McClintock, J. E. 2009, ApJ, 691, 847 
Li, L.-X., Zimmerman, E. R., Narayan, R., & McClintock, J. E. 

2005, ApJS, 157, 335 
McClintock, J. E., & Remillard, R. A. 2006, Black hole binaries, 

ed. M. Lewin, W. H. G. & van der Klis, 157-213 
Miller, J. M., Fabbiano, G., Miller, M. C, & Fabian, A. C. 2003, 

ApJ, 585, L37 

Miller, J. M., Fabian, A. C, & Miller, M. C. 2004a, ApJ, 614, 
L117 

— . 2004b, ApJ, 607, 931 

Mitsuda, K., et al. 1984, PAS J, 36, 741 

Narayan, R., & McClintock, J. E. 2005, ApJ, 623, 1017 

Novikov, I. D., & Thorne, K. S. 1973, in Black Holes, ed. C. De 

Witt & B. DeWitt (New York: Gordon & Breach), 343 
Penna, R. F., McKinney, J. C, Narayan, R., Tchekhovskoy, A., 

Shafee, R., & McClintock, J. E. 2010, MNRAS, 408, 752 
Remillard, R. A., & McClintock, J. E. 2006, ARA&A, 44, 49 
Sadowski, A., Abramowicz, M., Bursa, M., Kluzniak, W., Lasota, 

J., & Rozanska, A. 2010, ArXiv e-prints 
Shafee, R., McClintock, J. E., Narayan, R., Davis, S. W., Li, 

L.-X., & Remillard, R. A. 2006, ApJ, 636, L113 
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 
Socrates, A., & Davis, S. W. 2006, ApJ, 651, 1049 
Soria, R., & Kuncic, Z. 2008, in American Institute of Physics 

Conference Series, Vol. 1053, American Institute of Physics 

Conference Series, ed. S. K. Chakrabarti & A. S. Majumdar, 

103-110 

Soria, R., Risaliti, G., Elvis, M., Fabbiano, G., Bianchi, S., & 

Kuncic, Z. 2009, ApJ, 695, 1614 
Soria, R., Zampieri, L., Zane, S., & Wu, K. 2010, MNRAS, 1537 
Steiner, J. F., McClintock, J. E., Remillard, R. A., Gou, L., 

Yamada, S., & Narayan, R. 2010, ApJ, 718, L117 
Steiner, J. F., Narayan, R., McClintock, J. E., & Ebisawa, K. 

2009, PASP, 121, 1279 
Tsujimoto, M., et al. 2011, A&A, 525, A25-I- 

Webb, N. A., Barret, D., Godet, O., Servillat, M., Farrell, S. A., 

& Gates, S. R. 2010, ApJ, 712, L107 
Wiersema, K., Farrell, S. A., Webb, N. A., Servillat, M., 

Maccarone, T. J., Barret, D., & Godet, O. 2010, ApJ, 721, L102 
Zhang, S. N., Cui, W., & Chen, W. 1997, ApJ, 482, L155-I- 



